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Intermediate spin, which occurs in the theory of anyons, 
can also be exhibited by mesoscopic magnetic particles. The 
necessary broken time reversal symmetry is due to a suit- 
ably directed magnetic field. As a function of this field a 
system passes periodically through points which correspond 
to whole or half-integer spin. Intermediate spin is defined by 
fields which lie between these points. Since the tunnel split- 
ting in the ground doublet vanishes for half-integer spin, this 
splitting becomes periodic. The doubly periodic oscillations 
observed in the magnetic molecular cluster Fes represent the 
first unequivocal observation of this phenomena. Here the 
whole or half-integer nature of the system, i.e., the parity, is 
periodic for a field which is along either the easy or hard axis 
or a suitable combination of the two. A detailed theory is 
presented. 



The search for quantum effects in mesoscopic systems 
has become quite fashionable. In particular, relatively re- 
cently advances in several technologies have encouraged 
the investigation of the cross-over between the classical 
and quantum regimes in magnetic molecular clusters. A 
particularly interesting quantum effect which can occur 
in such large spin clusters has been emphasized by Loss 
et al. jl| and von Delft and Henley [0. For half-integer 
spin the amplitude for tunneling with a complete reversal 
of the magnetization of a small ferromagnet is zero due 
to interference effects arising from the so called topologi- 
cal term in the effective action. This parity effect implies 
that, while a whole integer spin system can have a tun- 
nel split ground state doublet, the half- integer analogue 
will not. Garg has shown, for a small easy axis ferro- 
magnet with a field perpendicular to the direction of the 
magnetization, the same tunneling amplitude oscillates 
as a function of field. 

The author has shown |0 that mesoscopic magnets, in 
a suitably directed magnetic field, can exhibit interme- 
diate spin. Intermediate spin values occur naturally 0] 
in the theory of anyons and are usually associated with 
the SO (2) group algebra which contains the single oper- 
ator S z with the spectrum S z — n — a/2, where n an 
integer, and where a is the usual statistical parameter. 
Evidently, anyons imply broken time reversal symmetry 
and, in the present context of a real physical magnet, 
the presence of a magnetic field. The necessary reduc- 
tion in symmetry from SO(3) (or SU(2)) to SO(2) (U(l)) 
will occur for large spin and an easy plane magnet. For 
such a magnet, it has been shown M the parity changes 



smoothly and periodically as a function of the magnetic 
field. The absence of a tunnel splitting for half-integer 
points directly implies oscillations in the tunnel splitting, 
i.e., the oscillations of Garg are a trivial consequence of 
intermediate spin. 

Very recently Wernsdorfer and Sessoli have experi- 
mentally observed double oscillations in the tunnel split- 
ting for Fes which is a molecular cluster with a net spin 
5 = 10. Oscillations which occur with a field along the 
hard axis are interpreted in terms of intermediate parity 
(i.e., spin) while those which are induced by a field along 
the hard axis are discussed in terms of the Garg [0 inter- 
ference effect. The principle purpose of this contribution 
is to show that both effects find a natural interpreta- 
tion and a relatively simple quantitative description in 
terms of the author's [Q approach to intermediate spin 
in magnets. It will be shown, in the context of Fes, that 
periodic changes in parity can occur even for an easy (y-) 
axis magnet when the full SO(3) or SU(2) algebra is im- 
portant. However, it must be that the symmetry in the 
plane perpendicular (x — y-) plane is broken, so that it 
remains the case that there is an easy plane for tunneling 
between the classically denned minima. 

It has been shown by numerical simulation [0 that Fes 
is adequately described by the Hamiltonian: 



H = -DS Z 2 + E (S x 2 - S y 



(1) 



where D > E so that the z-axis is easy while the x-axis is 
hard. The term involving S y 2 might be eliminated using 
Sy 2 = 5(5+1) — (S x 2 + S z 2 ) so that, to within a constant, 
H = -(D - E)S 2 + 2ES X 2 + S-h where h = gjJt B H is 
the field in energy units. 

Formally this Hamiltonian has been considered in some 
detail in and by Preda and Barnes using the auxil- 
iary particle method [0. This approach will be used for 
the detailed calculations. 

First consider a longitudinal field hi, i.e., one along the 
z-axis, and to be specific physical whole integer spin value 
5. The Hamiltonian only connects states \S Z — 5 — n >; 
n an integer, for which An is a multiple of two. It follows 
that when the levels S z = —5 and S z = 5 — n cross 
there can only be a tunnel splitting if n is even. The 
crossings at fields h = sD, with s even correspond to 
whole integer physics while those with s odd behave as if 
the spin was half- integer, i.e., a longitudinal applied field 
causes a periodic change of parity. 

At a more formal level, following ref. , it is observed 
that the longitudinal field hi can be eliminated by re- 
defining the spin operator to be 
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(2) 



which clearly has the intermediate spin spectrum S z = 
to — §; to an integer. This transformation does not 
change [S' Z ,S' ± ] = ±S' ± fcwi is not consistent with 
[S f+ ,S f_ ] = 25*2. The argument that the problem maps 
to itself is more subtle than in ref. jl. 

The periodic behavior as a function of a transverse field 
can be understood in terms of the fixed point E = D. 
Here, to within a constant, Ti reduces to +2ES X 2 — S x h t 
where h t is the transverse field in energy units. This 
defines another intermediate spin, or level crossing prob- 
lem. For an integer S, the ht = ground state is non- 
degenerate, while it is doubly degenerate corresponding 
to a half-integer point when ht = 2E. The "splitting" is 
periodic with period Ah t — AE. This fixed point behav- 
ior explains qualitatively the tunnel splitting oscillations 
as a function of /it, however a calculation of the general 
result for the period, the tunnel splitting, and the com- 
bined effects of longitudinal and transverse fields requires 
a detailed formalism and some tedious calculations. 

In outline the calculation goes as follows. The axis 
of quantization is put, x <-> z, along the direction of 
the transverse field direction. A basis \S Z >= \n > is 
chosen, and an auxiliary particle @,^|, a fermion f n , is 
associated with each state via the mapping \n >— > /t| > 
where | > is a non-physical vacuum without any auxiliary 
particles. Defined is a bi-quadratic version of an operator 
6 via: 6 -> J2n,n> fl < n\6\n' > /„,. The constraint 
Q = E„ = E„ fnfn = 1 - This procedure gives 



W = E„ - h t n) 

\D[{M^) 2 + (Mr 1 ) 2 ])/]/ n 



~ jDEn C +1 K + + l(4 2 /n + tifn+2) 
+ \h t En + fUn+ll 



(3) 



where the M™ +1 = [S'(S'+l)-n(n+l)] 1 / 2 are the matrix 
elements of S . Without a longitudinal magnetic field 
he, this is two isolated tight binding chains of spinless 
fermions /t. A longitudinal field hi provides a coupling 
between these two chains. The constraint Q = 1 implies 
this is a single particle problem. Notice that the "site" 
indices are either whole of half-integer following the par- 
ity of the spin. 

Schrodinger's equation 

(e — (2En 2 — h t n))a n = 

+ \D[MZ+ l NQXla n+2 + MZ+Hl^a,^ 

+ [(Af«+ 1 ) 2 + (Mr 1 )V 

- i^(M™ +1 a„+i + M£-V-i), (4) 

involves finite differences, where the wave-function ^ = 
En a ™/nl >■ The difference in the site energies is ~ 2E 



near the center n ~ of the chains while the "hopping 
term" is ~ S 2 D. For the problem of interest D is positive 
and the ground state wave function does not change sign. 

The ground state is located near the center of the 
chains and for most purposes it will suffice to replace 
the M™ +1 by a constant matrix element S whence the 
problem reduced to that of a (discrete) harmonic oscilla- 
tor. It might be expected that continuum approximation 
is valid since S 2 D 3> 2E. If, he — 0, and if the passage 
is made from the discrete to a continuum problem the 
two chains are are identical and there is no tunnel split- 
ting. In this rather unusual representation of a tunneling 
problem it is the exponentially small difference in energy 
between the ground state on the even site chain (the ab- 
solute ground state) and that on the odd site chain which 
leads to the "tunnel" splitting. 

However, in the absence of fields and if E = this 
equation must reproduce the eigenenergies —Dm 2 where 
to = —S — (S — 1) . . . S — 1, S is an integer quantum 
number. In particular the ground state is exactly doubly 
degenerate with to = ±5 and again with one such state 
on each chain. In order to obtain correctly this limit use 
is made of the exact result for E = H = 0, namely 
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(2S)\ 



{S-n)\(S + n)\ 



1/2 



(5) 



where the two solutions correspond to S z = ±S. For 
large spin S this approximates well to 



2 s 



(6) 



i.e., ignoring the factor e^ l7r the wave-function is that of 
a simple harmonic oscillator fairly well localized relative 
to the limits at n — ±S. It can be made to lie on one 
or the other chains by taking the sum or difference of 
these solutions. The system does not "see" that n is lim- 
ited. Adding the term 2En 2 causes the wave function to 
become slightly more localized. Since this additional po- 
tential is harmonic the wave function remains harmonic. 

Consider the problem with hg = but with both E 
and h t finite. The wave function is taken to have the 
form 



f(ri)a; 



±s 



(7) 



Substituting this into Schrodinger's equation gives at a 
sufficient level of approximation: 



e/(n) = (2En 2 - h t n)f{n) 
+ \S 2 D[f(n + 2) + f(n-2)-2f(n)}. 



(8) 



At this point the intermediate spin transformation is 
made, i.e., the field h t is absorbed into the term 2En 2 
by the shift n — * n — (ht/iE). There is an appar- 
ent displacement of S = (ht/AE). The small change 



2 



in energy is unimportant and will be dropped. How- 
ever the real displacement d is different. With h t = 
to a good approximation, f(n) ~ e -(2E/D) 1/2 n 2 /2S wn jj e 
for finite h u f(n) ~ e ^E/D)^ {n ^ ht/iE)r /2S sq ^ 

a n ~ e -™ 2 /2S e -(2£/D) 1 / 2 („-(h t /4£)) 2 /25 which ig approx _ 

imately e ( n -(' I t/ 2 ( 2BD ) 1/2 )) 2 / 2S ] i.e., the real displacement 
of the net solution is 

d=h t /2{2ED) 1/2 , (9) 

which corresponds to a statistical parameter a = 
ht/ (2ED) 1 / 2 . The point to be made is that to a good 
approximation the solution for finite h t is the same as 
for ht = except for a displacement of the site location 
from integer (or half-integer values for half-integer S) by 
d. 

A Fourier transform, is made in order to examine 
the approach to the continuum limit and to generate 
a tunnel-barrier problem: j(y) = -^== J dn e my a(n) , 
whence Schrodinger's equation (he = ht = 0) is: 

( 2E ^2)f(y) = "^[cos2y - l}f(y), (10) 

which is Mathieu's equation ||. The potential is periodic 
with period Ay — it and the solutions form bands. For 
a given energy, a solution might be characterized by a 
wave- vector k and Floquet's (Bloch's) theorem (9) implies 
that solutions are of the form fk{y) — e lky Uk(y) where 

Uk{y) = Mfc(p + 7T). 

Since the fc-space is reciprocal to the y-lattice it is 
roughly equivalent to the real x-space and since this lat- 
ter is discrete, making an inverse Fourier transform se- 
lects a particular wave vector k. Also a shift in fc-space 
is fully equivalent to a shift in real space. Because of the 
intermediate spin shift d, the sites are precisely shifted 
from the whole or half integers. Specifically for a whole 
integer spin, the solution which is finite on the even site 
chain implies k = d. Similarly for the odd site chain 
k = d — 1. The only difference for a half-integer spin 
is an addition shift of 1/2. Without a intermediate spin 
(i.e., statistical shift) this makes the two k- values equiv- 
alent at k = ±1/2 and directly implies a zero splitting in 
the absence of fields. 

In the y-space, around y = 0, f(y) = 
(l/ v / /V^)e- s/2 / 2/32 with /3 2 = S^TEjD, and the nomi- 
nal ground state energy is (wo/2); ujo = 2Syj2ED. The 
band energies are: tk — (^o/2) + {w/2) cos7rfc where the 
width H w = 8^2/TTijJoSy 2 e~ s - f and where the action 
S f = 2S y / {D/2E). The result for the "tunnel splitting", 
i.e., the difference in energy between the ground and first 
excited states is 

5E = ^^Lo S) ,2 e- s f |cos(7rd)| ; 

S f = 2S^D/2E). (11) 



This result for the tunnel splitting is very similar to that 
for an easy plane magnet [00]. The most significant 
difference is that the period for the oscillation, with a 
field along the hard axis is 2(2ED) 1 / 2 rather than, in 
the present notation, 4E. For E small this result for the 
period agrees with the result of Garg 2(2E(D + E)) 1 / 2 
quoted by Wernsdorfer and Sessoli Q . 

The last step is to account for the longitudinal field 
hg. There is a level crossing for the values hi = sD; s 
an integer. Since they reflect the essential physics only 
results for these crossing points will be presented. Con- 
sider a whole integer physical spin S, following the earlier 
discussion there are two types of crossing corresponding 
to the value of s. For s even, including s — 0, the system 
behaves as if it is whole integer while for s odd it has the 
characteristics associated with a half-integer total spin. 
The proof goes as follows: The zero order wave functions 
for the states which cross are, e.g., azf and af _s which 
are both exact solutions. FromEqn. (||): a~ s = (— l)"af 
where is even in n, while a^~ s is even or odd as s is 
even or odd. Given there is tunneling, the appropriate 
wave functions are of the form, for large enough S: 

a ± , n = f(n)[(-l) n a^±at s }. (12) 

If this function is substituted in Eqn. (^) the terms in- 
volving hi cancel to O(S'), so that the resulting difference 
equation and solution for f(n) are unchanged. However, 
the function a± j?l is quite different depending on the even 
or odd character of s. With f(n) = 1, because of the 
term (—1)™, a± jT1 jumps wildly between sites but defines 
a smooth function on either the even, „, or odd , 
a ° ± n , sites. When s is odd it is almost the case that 
a e ± n — —a° ± _ n — — a^j- n . In fact, if the smooth function 
which passes through the lattice points is displaced by 
d = Sn = 1/2, then these relationships are exact. As dis- 
cussed above, it is just such a displacement with converts 
the whole into the half integer problem and demonstrates 
that the s odd crossings are equivalent to a physical half- 
integer spin S — s/2. When s is even the functions „ 
and dj_ n are very different, however (to within a sign) 
a§_ n = a^- n so that difference between the two solutions 
a± tn is a displacement by one lattice site or d = 1 which 
is just the generalization of the same result for s = 
and which results in a the maximal tunnel splitting with 
h t = 0. The tunnel splitting reflects the decrease in ef- 
fective spin in the sense that it increases strongly with 
•s. 
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